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Topological insulators embody a new state of matter characterized entirely by 
the topological invariants of the bulk electronic structure rather than any form 
of spontaneously broken symmetry jl-8|. Unlike the quantum Hall or quan- 
tum spin-Hall-like systems, the three dimensional topological insulators can host 
magnetism and superconductivity which has generated widespread research ac- 
tivity in condensed-matter and materials-physics communities. Thus there is an 
explosion of interest in understanding the rich interplay between topological and 
the broken-symmetry states (such as superconductivity) j9l-l2C)(|. greatly spurred 
by proposals that superconductivity introduced into certain band structures will 
host exotic quasiparticles which are of great interest in quantum information sci- 
ence. The observations of superconductivity in doped Bi2Se3 (Cua;Bi2Se3) and 
doped Bi2Te3 (Pdx.-Bi2Te3 Tc ~ 5K) have raised many intriguing questions about 
the spin-orbit physics of these ternary complexes while any rigorous theory of 
superconductivity remains elusive [211, l22l | . Here we present key measurements 
of electron dynamics in systematically tunable normal state of Cua;Bi2Se3 (x=0 
to 12%) gaining insights into its spin-orbit behavior and the topological nature of 
the surface where superconductivity takes place at low temperatures. Our data 
reveal that superconductivity occurs (in samples compositions) with electrons 
in a bulk relativistic kinematic regime and we identify that an unconventional 
doping mechanism causes the topological surface character of the undoped com- 
pound to be preserved at the Fermi level of the superconducting compound, 
where Cooper pairing occurs at low temperatures. These experimental obser- 
vations provide important clues for developing a theory of superconductivity in 
highly topical Bi-based 3D topological insulators. 

Bismuth selenide has been experimentally discovered by angle-resolved photoemission 
spectroscopy (ARPES) to be a topological insulator with a large bulk bandgap 
Spin- resolved photoemission studies reveal that surface electrons in Bi2Se3 form a Dirac 
cone spanning the bulk insulating gap, composed of spin-momentum-locked helical states 
. Recent theory suggests that just as topological insulators can be identified from band 
structure alone, the topological properties of spin-orbit superconductors can be determined 
from the connectivity of bulk and surface band structure 10l-ll3l|. Superconducting pairing 



induced on a two dimensional topological surface state via proximity to a bulk supercon- 
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ductor was shown to be inherently unconventional, with the potential to host non-Abelian 



vortices 
ing 



20| and defying traditional classification as a singlet or triplet form of Cooper pair- 



12| . To ensure that a material is capable of hosting two dimensional vortices in the 



superconducting state, the spin-polarized surface states are must be non-degenerate with 
the bulk bands where Cooper pairing takes place. It has been found that 10% copper is 
needed to bring about superconductivity in bulk Cua;Bi2Se3, with a relatively large tran- 



sition temperature of Tc=3.8° K confirmed by the observation of the Meissner effect [21|. 
All known measured band structures as well as band calculations of undoped Bi2Se3 suggest 
that the topological surface states cannot remain non-degenerate at the Fermi level under 
the carrier doping that is required to induce superconductivity in Cu2;Bi2Se3. Therefore it 
was expected that the topological effects on the surface of Cui,.Bi2Se3 would be washed out 
at the superconducting compositions (comparisons in online Supplementary Information). 

In clear contrast, as we show here a direct band structure measurement of supercon- 
ducting Cuo.i2Bi2Se3 reveals that the surface states are distinctly isolated from the bulk 
conduction bands where bulk superconductivity takes place (see Fig. 3). Motivated by this 
puzzling observation we have systematically investigated the doping process by examining 
several stoichiometric crystals with copper added (Cu2.Bi2Se3, x=0, 0.01, 0.05, 0.12) and 
contrasting these doping levels with a case in which copper is substituted for bismuth as 
Cuo.iBii gScs. These results, summarized in Fig. Ic, suggest that copper atoms can add 
holes or electrons depending on their net preference for occupying different kinds of site in 
the Bi2Se3 lattice. Copper atoms intercalated between van der Waals-like bonded selenium 
planes {Cuint) in the crystal are thought to be single electron donors, while substitutional 
defects in which copper replaces bismuth in the lattice (Cu^i) contribute two holes to the 
system which is similar to the behavior reported in the impurity atom study on Bi2Se3 
matrixj2^. The qualitative result of copper addition is a gradual enlargement of the ob- 
served Fermi surface from electron doping and, further surprisingly, a strong reduction in 
surface electron velocities through reshaping of the surface state conduction band. The slope 
of the surface state band represents particle velocity, and is globally reduced in amplitude 
by 30% above the Dirac point when copper doping of x=0.05 is added to the stoichiometric 
compound (see Fig. 2c). Assuming that the surface state Dirac point energy is fixed relative 
to bulk bands, reducing the slope of the surface state increases separation between the bulk 
and surface conduction bands, a characteristic that increases stability of the topological 
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surface state and is of great importance for stabilizing the topological nature of supercon- 
ductivity. This effect acts together with the nonlinear doping mechanism reported here to 
preserve the two dimensional topological surface state at superconducting doping. Therefore 
our results suggest that the crystal surface of Cux.Bi2Se3 can host two dimensional vortices, 
which is a prerequisite for observing non-Abelian statistics on the surface. 

Addition of copper causes the surface state kinetics to become strongly hexagonally 
anisotropic, (Fig. 2b) a deformation that makes the surface susceptible to spin-fluctuation 



or magnetic instabilities [17||. This is due to the fact that the topological surface states in 
Bi2Se3 are spin polarized [6|. Carrier density in the surface state is much greater than den- 
sity in the bulk, as estimated by the Luttinger count applied to the photoemission (ARPES) 
data (Fig. Ic), suggesting that the surface carries a screened negative charge within the 
normal state of the superconductor. We also observed that attempting to force the creation 
of Cusi replacement defects by adding less bismuth results in very weak hole doping for 
Cuo.1Bi1.9Se3 (Fig. Ic), raising the bulk conduction band entirely above the Fermi level. 
Upon doping into the superconducting regime (x=0.12), the Fermi energy is found to be 
raised to 0.25 eV above the bulk conduction band minimum, placing the Fermi level in the 
linear (relativistic) regime of bulk band dispersion. Contrasting the bulk band lineshape 
with classical and relativistic Dirac-like energy dispersions with the same mass (M=0.155 
Me) reveals that electron kinetics begin to enter the linear relativistic regime within ~0.1 
eV of the Fermi level (Fig. 3f). A slight bend in the dispersion centered near 90 meV 
binding energy (Fig. 3f inset) may suggest electron-boson interactions in the system, also 
ubiquitously observed in other superconducting materials. Therefore, Cooper pairing in the 
superconducting state takes place in this relativistic regime where the chemical potential 
lies. A three dimensional massive Dirac-like dispersion in topological spin-orbit materials 
is expected as a direct result of the band inversion mechanism that causes the topological 
insulator state, and has also been experimentally observed in other topological insulator 
Bii_2^Sba; alloys [4]. In Cu2;Bi2Se3, calculations predict that band inversion occurs at the 
F-point {kx=ky=kz=0) in the center of the bulk conduction band leading to a spin-orbit 
induced Dirac-like bulk band in clear qualitative agreement with our experiments. 

Several features of this unusual spin-orbit band structure provide critical insights into 
characteristics of the superconducting wavefunction. Bulk Fermi momenta of 0.110±3 
and 0.106it3 are observed along the F-M and T-K directions respectively. Varying 
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incident energy to observe dispersion along the z axis (F-Z direction) reveals a Fermi mo- 
mentum of 0.12±1 A~^, suggesting that the bulk electron kinetics are three dimensionally 
isotropic. Carefully tracing the band (Fig. 3e,f) reveals a Fermi velocity of 3.5±2 eVxA 
along F-M and 4.1 ±2 eVxA along F-i^, estimated within 50 meV of the Fermi level. The 
gap between bulk valence and conduction bands appears to be unchanged upon copper dop- 
ing. In a material specific band-structure based modeling of superconductivity, considering 
a competition between inter-site and intra-site orbital coupling suggests that the supercon- 
ducting pairing symmetry in this system is partially determined by a ratio between the Dirac 
rest mass and doped chemical potential ("/x") tl2|], when the conduction band is fitted to a 
standard relativistic Hamiltonian: 



E^Oz) = MTo + v*ciKT, + kyT2 + fc.Fg), (1) 

where M* is the rest mass, vc is the velocity of the 3D Dirac band, and the Fj's (i = 
0; ... ; 3) are the 4x4 Dirac Gamma matrices. For this discussion we assume that the 
Vc is isotropic, although our data suggest that it may be smaller in the z-axis direction. 
The chemical potential of this Hamiltonian is equal to the sum of the rest mass and the 
binding energy of the valence band minimum, which can be experimentally evaluated from 
our photoemission data. In an ideal Dirac band structure, the rest mass will be equal to 
half the band gap, or 150 meV (giving M*=.15eV, /i=0.15 eV + 0.25 eV, and ^ ~ |). Due 
to asymmetries between conduction and valence band dispersion in the real experimental 
system, a larger effective rest mass is used to fit the curvature of the bottom of the conduction 
band in Fig. 3f, consistent with a ratio of ^ =|. These ratios are small enough to allow the 
rare theoretical possibility that the bulk superconducting wavefunction may have an odd (-1) 



parity symmetry value [12|]. A dispersion anomaly that we have noted in the inset to Fig. 3f 
also effects low energy electron kinetics, making the bands more linear locally at the Fermi 
level. Dispersion anomalies, often referred to as "kinks" in photoemission literature, are a 



common property of band structure in real superconducting materials 



25| . The presence of 



the kink depresses the value of — in the favorable direction for realizing odd parity pairing 



proposed in Ref. |l2 |. 



The observed spin-orbit band structure kinematics fixes several key parameters related 
to the superconductivity. The superconducting coherence length can be estimated from 
the average Fermi velocity we observe here to be about 2000A (,^o~ 0.2 x Tivf/KbTq = 



0.2 X 3.8eVA/{KB x 3.8°-?^) = 2000A), assuming minimal scattering and a superconducting 
gap related to Tq by the BCS formula. This value is about 1-2 orders of magnitu de g reater 
than that seen in the cuprates (^o~ 100 — 200A) or cobaltates (^o~ 200A) 25|, l26[. The 



phase ordering temperature scale, also estimated from our data, is approximately Tg = 
60,000°K, four orders of magnitude larger than the superconducting critical temperature of 
3.8°K (Te= ^0 x h^nj2m*= 2000A x x 10'^'^cm-y{2 x .155Me) = 60,000°K). A large 
coherence length and high phase ordering temperature suggest that key properties of the 
pairing state can be described within a mean field picture. Based on this, we expect the 
superconducting gap to be about ~0.6 meV (3.5xKbTc/2=0.6 meV). This is much smaller 
than the state-of-the-art resolution of ARPES. 

The fact that superconductivity occurs in the relativistically linear bulk bands in 
Cu2:Bi2Se3 as seen in our data allows unusual combinations of spin and orbital mixing in 
the bulk superconducting state, making pairing take an unusual form; more remarkably 
however, probing near the Fermi level, our data show that the spin-textured surface and 
doubly degenerate bulk bands are well separated by 0.04 in momentum and an energy 
spacing of about A£;=130 meV, establishing that bulk superconducting pairing occurs in 
the presence of a non-degenerate, spin polarized two dimensional topological surface state, 
and that theoretical constraints ensuring unconventional superconductivity induced in the 
;opological surface band structure will define key properties of the pairing wavefunction 
3], Q. The remarkable preservation of the surface state occurs due to nonlinear doping 
and unusual distortions of the surface state that we have discovered here (Fig. 2c). This 
experimental observation is in contrast to first principles predictions, which suggest that 
the surface state will become degenerate and thus lose the protected spin-texture at the 
Fermi level due to intersection with the bulk bands (see Fig. 4b, Section SI III in online 
Supplementary Information). 

Topological surface states are completely determined by the bulk band-structure topology 
[3I, ensuring that superconducting symmetry breaking in the bulk bands must necessarily 
effect the surface spectrum. Three dimensional topological insulators doped into a supercon- 
ducting phase are classified by an integer invariant ("n"), rather than a Z2 invariant (parity 
invariant) [l2|. The specific band structure and strong spin orbit coupling of Cu2;Bi2Se3 
suggest that the topological invariant "n" may be non-zero, allowing the possibility of topo- 
logical superconductivity in Cua;Bi2Se3. This is further supported by the fact that the ratio 
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^ for Cu2;Bi2Se3 is well into the topological regime [12| as seen in our data. The exact 
nature of the superconducting order-parameter will depend on the parity of the pairing 
wavefunction, which may be even or odd under inversion since Bi2Se3 is a centrosymmetric 
crystal. Of the two remaining possibilities, even parity will result in a fully gapped band 
structure, and odd parity will generate new surface states resembling those shown in Fig. 



4d 1^. In either scenario, the preserved topological spin polarized surface state will allow 



the superconducting wavefunction to host Majorana fermions [20|, which might be manip- 
ulated adiabatically if the surface is fully gapped (even parity. Case 1 in Fig. 4). In the 
second, odd parity case, the material will be the first known realization of a state of matter 



specifically termed "topological superconductivity" 12 



18 



While the exact determi- 



nation of order-parameter phase will require low temperature phase sensitive measurements 
on ultra-high purity large single crystals, the unusual doping evolution towards achieving 
superconductivity discovered here provides critical clues for understanding pairing in doped 
topological insulators, as well as constituting the key ingredients for developing a general 
theory of these novel Bi-based superconductors. 
Methods summary: 

Experimental observation of surface and bulk separation and their microscopic electron 
kinetics constitutes overcoming several significant technical challenges, (see online Supple- 
mentary Information) 

Angle resolved photoemission spectroscopy (ARPES) measurements were performed at 
the Advanced Light Source beamlines 10 and 12 using 35.5-48 eV photons and Stanford 
Synchrotron Radiation Laboratory (7-22eV photons) with better than 15 meV energy reso- 
lution and overall angular resolution better than 1% of the Brillouin zone (BZ). Samples were 
cleaved and measured at 15°K, in a vacuum maintained below 8x10^^^ Torr. Momentum 
along the z axis is determined using an inner potential of 9.5 eV, consistent with previous 
photoemission investigations of undoped Bi2Se3 Large single crystals of Cua:Bi2Se3 were 
grown using methods described in the supplementary information. Surface and bulk state 
band calculations were performed for comparison with the experimental data, using the 
LAPW method implemented in the WIEN2K package {2^. Details of the calculation are 
identical to those described in ref. Q|. 
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FIG. 1: Superconductivity in Cua.Bi2Se3 crystals: a, A low-energy electron diffraction im- 
age taken at 200eV electron energy provides evidence for a well ordered surface with no sign of 
superstructure modulation, b, Resistivity and magnetic susceptibility measurements for samples 
used in this study, (detailed materials characterization is reported in Ref. 2l[) Samples exhibit 



a superconducting transition at 3.8°K at optimal copper doping (x=0.12). c, The number of 
charge carriers is calculated from the Luttinger count (^^f^ff i x2 for the doubly degenerate bulk 
band). LDA predictions show the carrier density obtained by aligning LDA band structure with 
the experimentally determined binding energy of the Dirac point. 
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FIG. 2: Surface electron kinematics of Cua;Bi2Se3: a, The hexagonal surface Brihouin zone of 
Cua;Bi2Se3 is shown above a diagram of the three dimensional bulk Brillouin zone, b, Symmetrized 
surface state Fermi surfaces are displayed with (right) a GGA prediction based on the experimental 
Fermi level, c, Photoemission measurements are shown at high photon energy (E>20 eV) for 
non-superconducting Cua;Bi2Se3, demonstrating a reduced surface state dispersion after copper is 
added, d, When the Bi-Se plane spacing at the surface of a 12-layer slab is increased by 0.2 A 
(1.64 to 1.84 A), dispersion in the upper Dirac cone increases by 16%. 
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FIG. 3: Band topology at the superconducting composition: a-b, Momentum dependence 
of the bulk and surface conduction bands in superconducting Cuo.i2Bi2Se3 is measured with low 
energy (9.75eV) photons for enhanced bulk sensitivity, c, Dispersion along the z-axis is examined 
by varying incident photon energy, d, Bulk and surface bands remain separate at intermediate 
kz values, e-f, Energy of the bulk electrons is compared with Dirac (vc=6 eV-^) and classical 
(parabolic) fits with a mass of 0.155Me. An inset shows self-energy with respect to the Dirac fit. 
g, The surface electronic structure presents a non-trivial topological setting for superconductivity 
because (green) surface and (blue) bulk bands do not overlap. 
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FIG. 4: Symmetry breaking on the topological surface: a, Bulk and surface electrons are 
non-degenerate at the Fermi level, b, A phase diagram compares (expt) the measured supercon- 
ducting topology with preliminary expectations based on cases in which (A) each Cu atom donates 
one doped electron, (B) the experimental chemical potential is applied to GGA band structure, 
and (C) the experimental chemical potential is applied to the band structure of undoped Bi2Se3. 
Electronic states expected below Tc are drawn for (c,e) even parity superconductivity and (d,f) 
an example of odd parity "topological superconductivity" . Panels e-f show states within 5 meV 
of the Fermi level. 
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